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1. Introduction 

For a few years we have had algorithms that allow simulation at small quark masses m in 
large volumes, such that M^^L > 4. As a consequence of this progress, and since volume effects 
are exponentially suppressed in Mj^L, they are presently kept at the level of 1%. On the other 
hand, discretization effects obey power law scaling, which in the case of 0{a) improved Wilson 
fermions amount to 0{a^) effects. Most importantly, for all kinds of Wilson fermions, the dominant 
discretization effects are related to the breaking of chiral symmetry. As a consequence, they become 
more relevant as we try to reach lower quark masses at fixed a. It is then clear that the understanding 
and removal of discretization effects is essential. 

Recent studies [1] however have shown a severe critical slowing down which makes it unreal- 
istic to simulate lattices with a < 0.05 fm. While the critical exponent can be reduced from ~ 5 to 
2 by using open boundary conditions [2], to reduce the lattice spacing remains challenging and an 
understanding of the discretization effects is necessary to establish trustworthy results. 

One framework to study the discretization effects is provided by Wilson ChPT (WChPT) [3], 
the chiral effective theory of Wilson LQCD. In this effective theory the cutoff effects are parametrized 
at a given order in the chiral expansion by a finite number of Low Energy Constants (LECs), among 
which we find C2. In terms of the physical and unphysical LECs, Wilson ChPT predicts the scaUng 
of pion observables with respect to m and a. In this context, for one class of observables, C2 plays 
an important role as it appears in the LO lagrangian. 

In this talk we present our results and strategy to put at least an upper limit on \c2\ and even- 
tually establish its sign. We also discuss the relevance that the values found would have on the 
m-dependence of some physical observables. 

2. Wilson ChPT 

The first step to build WChPT is to collect all the operators in the continuum Symanzik effec- 
tive theory that are compatible with the symmetries of Wilson LQCD up to a certain dimension (i.e. 
power in a). These then have to be matched to the chiral effective theory; they must be rewritten 
in terms of pion fields. Once a power counting for a relative to m and p has been established, the 
spurion technique which was already used to introduce chiral symmetry breaking effects due to m 
can be extended to treat cutoff effects. 

Our aim is to demonstrate that cutoff effects are not strong enough to distort the physical chiral 
scaling of the observables in which we are interested. As a consequence, we measure on our coarser 
lattices with a = 0.075fm and a = 0.065fm and assume a conservative power counting for a (the 
so called LCE regime) which assumes: 

m ~ a^AgcD a 2± 0.07 fm ^ a^A^^D ~ 10 MeV . (2. 1) 

This counting is expected to be the proper one for quark masses around lOMeV. However large 
or small coefficients in the expansion could shift this limit up or down so a practical test is nec- 
essary. We also assume (which is true in our case) that 0{a) improvement has been carried out 
nonperturbatively. 
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Ml = Ml + 0(m4) mI = -^- 2c2fl2 (2.3) 



With the power counting in eq. (2. 1), the LO chiral Lagrangian is: 

= — Tr [d^Ud^U^] - — Tr [U + U'] + ^y^Tr [U + [7^]^ (2.2) 
where C2 is a new unphysical LEC. One can show that the pion mass is corrected at LO: 

and that a new phase pattern emerges [4], whose properties are controlled by the sign of C2 [3]. For 
Wilson fermions, if C2 > it is possible to reach a phase where flavor and parity are spontaneously 
broken by decreasing the quark mass at fixed a. Since the transition is of second order, the pions 
ai^e massless on the boundary. Conversely, if C2 is negative, the transition is first order, flavor and 
parity remain unbroken, and there are no massless pions at finite a. 

We now possess a tool to predict the scaling on m and a (up to a given order) in terms of the 
physical LECs of ChPT plus some new ones (including C2). It should be noted that m appearing 
in the Lagrangian (Eq. (2.2)) is just a bare parameter, which has no immediate relation to a well- 
defined quark (or pion) mass of the lattice theory. One example of such a quark mass is the PCAC 
mass, nipcAC- Before the WChPT formulas can be used, it is necessary to find a relation between 
m and nipcAC (or Mji) that allows to rewrite the observable of interest ff: 

G(m(mpcAc)) or d{M^) = G{m{M^) . (2.4) 



Ml=Ml 



M2 + 10c2a^ (Ml\ , J ; 



For example, one can obtain the pion mass as a function of Mq (or m through Eq. (2.3)) [5] ' : 

■0{a'), (2.5) 

'I ' 

but this formula cannot be used for extrapolations until a relation between Mq and a measurable 
mass is found. Unfortunately we know just the LO relation between m and mpcAC [6]: 

mpcAC=[m — \ [l + 0{a^ ,m)) . (2.6) 

After the adequate substitutions are performed, one finds that observables can be separated in two 
classes, according to how much they are affected by discretization effects. 

• Observables that are nonzero in the chiral limit, such as the pion decay constant. Cutoff 
effects appear here just as analytic contributions and at NLO, both in terms of M;^ and mpcAC- 



Ml Ml , -y , 



(2J) 



Observables that vanish in the chiral limit, such as the pion mass or scattering lengths. Cutoff 
effects appear at LO, they change the coefficient of the chiral logarithms in such a way that 
the expansion breaks down at low M^r- Here we show the S-wave scattering length for 1 = 2, 
Aq, as an example, that will also be useful in the following: 



2^ Mj 



3M2 + 12c2fl2 Ml 



32n^Fi 



lc2(r 
16nFl 



(2.8) 



'Here and in the following we do not write the analytic terms (polynomial in ) explicitly for simplicity, but they 
have also been calculated. 
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Figure 1: R(t) for the F7 ensemble. The blue curve represents the WChPT prediction. 



3. Determination of C2 

From Eqs. (2.6) and (2.3) it is apparent that in the relation Ml{mpcAc)^ C2 appears only at 
NLO. Conversely C2 enters at LO in the relation between the S-wave scattering length and Mj^, see 
Eq. (2.8). The energy eigenvalues of interacting particles are affected in finite volume by a term 
proportional to the scattering length according to the Liischer formula: 

E2„{p = 0) = 2M^ + AE AE = -^ + 0{L-^) (3.1) 

which opens up the possibility to measure C2 through volume effects. 
We consider the ratio of correlators: 

where 7r^(?) = Lx^^(^x). Apart from the QCD coupling and quark masses, R{t) needs no renor- 
malization and at LO in ChPT it depends only on Fji, Mji and C2. The result for R{t) on one of our 
lattices is shown in Fig. 1 . 

We have evaluated R{t) on the coarser CLS lattices listed in table 1 and in WChPT, in the 
LCE regime and finite volume. We evaluate all-to-all quark propagators stochastically, using 4 
fully time diluted noise sources per configuration. Using Fj^ and Mji determined from the two point 
correlation functions, we invert the ChPT formula to obtain an effective C2{t). As shown in Fig. 2 
we then average over a time interval [tminjmax] such that ?„„•„ > k/Mji, {T — tmax) > k/Mj[ for fixed 
A: = 0(1). In practice we have checked that the result does not change significantly for ^ G [2,3] 
and in the following we report the results obtained for ^ = 2.5. This definition converges to the 
LEC C2 in the chiral limit. We treat errors coming from autocorrelations in a conservative way, 
using the techniques described in [1]. To estimate the systematic effect coming from truncating the 
chiral expansion at LO we compare the results using in the chiral formulas both the observed pion 
decay constant F,^ and its value in the chiral limit Fq, see Fig. 2. 
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label 


a{fm) 


(MeV) 




A^cfg 


A4 


0.075 


380 


4.1 


100 


As 


0.075 


330 


4.7 


325 


E5 


0.065 


440 


4.7 


999 


Fe 


0.065 


310 


5 


300 


F7 


0.065 


270 


4.2 


600 



Table 1: Lattices on which R{t) was computed. Lattices used for physical measurements have lattice spac- 
ings as small as a = 0.05 fm. The lattice spacing was set from fx — 155 MeV, see contribution by M. 
Marinkovic at this conference. 
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Figure 2: The effective C2{t) with the observed decay constant Fj^ (left) and its value in the chiral limit 
(right). 

The results that we obtained are summarized in Fig. 3. Both determinations (using Fji and 
using Fq) seem to converge to a positive value, and all our points lie far from the dangerous region 
(the blue line) where the physical (due to M;^ / 0) and unphysical (due to a / 0) contributions to the 
S-wave 1 = 2 scattering length are equal. It is interesting to mention that the ETMC collaboration 
found a negative value [7], which is not in contradiction with our result, since their action was 
different. 

4. Conclusions and Outlook 

Our aim is to find whether or not we expect large cutoff effects in our simulations. Considering 
that we have performed tests on our coarser lattices and that we have assumed a pessimistic power 
counting, Fig. 3 is a good indication that cutoff effects are small in our observables and the physical 
scaling with (or nipcAc) is not significantly distorted for the masses that we are considering. 
Our procedure is cheap and simple, and will be repeated in the future for all ensembles. 

We have significant errors, coming mainly from the chiral truncation, that make it difficult to 
establish a conclusive value for C2, although our data at a = 0.065 fm point to a rather small and 
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F7, F6, E5 
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Figure 3: Results for C2 on lattices with a = 0.065 fm (left) and a = 0.075 fm (right). Thick lines corre- 
sponds to while thin lines corresponds to Fq. The dashed line indicates the physical point, the blue line 
corresponds to 2c2fl^ = mJ, which is where the physical contribution to the scattering length is of the same 
magnitude of the one coming from cutoff effects. 

possibly positive value. 

In the near future we will apply the same procedure to lattices with smaller masses than the 
ones so far considered. These additional points may help to establish the sign of C2- If indeed 
C2 > with our action, this would suggest adding a clover term to the action used by the ETM 
collaboration also yields C2 > [8]. This is important, because in that context Mj^± — Mjjo oc C2 and 
a negative C2 means large volume effects through Mj^a < M-^± . 
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